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Abstract 

We study the weakly-interacting Bose gas in both two and three dimensions 
using a variational approach. In particular we construct the thermodynamic 
potential of the gas to within ladder approximation and find by minimization 
an accurate mean-field description of the dilute Bose gas. Using spin-polarized 
atomic hydrogen as a specific example, we obtain an improved description of 
the Bose-Einstein condensed phase in three dimensions and a signature of a 
phase transition in two dimensions. The latter cannot be found by straight- 
forward application of perturbation theory around the ideal Bose gas. 
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I. INTRODUCTION 



The theoretical study of a three-dimensional dilute Bose gas in the regime where quan- 
tum degeneracy effects play a role has a long history, 0 and was in first instance important 
for a microscopic understanding of the A— transition in liquid 4 ife, although quantative 
predictions for this strongly-interacting system cannot be obtained by these means. The 
latter is however not true for magnetically trapped quantum gasesjl which explains why the 
phenomenon of Bose-Einstein condensation, accompanied by the occurrence of superfluidity, 
is still an active field of research today. Moreover, it turns out that there is no completely 
satisfying theory for these systems, !i since existing treatments, which are based on a Bogoli- 
ubov kind of approach, fail for instance to show the correct order of the phase transition and 
yield a first-order transition to the condensed phaseiJl instead of the second-order transition 
expected from the theory of critical phenomena. In this paper we will develop a variational 
approach and show for example that only a slight change in the Bogoliubov theory can 
resolve this issue. 

The two-dimensional dilute Bose gas has become subject of study only more recently 
when it was realized that this system is in the same universality class as the XY-model0 and 
should therefore exhibit a Kosterlitz-Thouless phase transition.! Physically this implies that 
there is no macroscopic occupation of one single quantum state, but there is nevertheless a 
transition to a superfluid state through the mechanism of binding of vortices with opposite 
vorticity. This mechanism causes the superfluid density n s to show a universal discontinuity 
such that n s A 2 = 4 at the critical temperature, i where A = {2’kTl 2 /rnksT) 1 / 2 is the thermal 
de Broglie wavelength and m is the mass of the particles in the system. The study of the two- 
dimensional Bose gas starting from a microscopic Hamiltonian was initiated by Popov.0 In 
his approach, which corresponds to a one-loop calculation and is based on the introduction 
of a so-called quasicondensate, the superfluid density exhibits no jump, but instead goes to 
zero in a continuous fashion. In addition, it has been showrffi that the range of validity 
of Popovs approach is by no means attainable in the experimental set-ups envisaged at 
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present. Therefore, we have recently extended this work using the ladder or T— matrix 
approximation which should be sufficiently accurate at low densities since it includes all two- 
body processes.0 In this manner we correctly found a jump in the superfluid density, but 
unfortunately, due to infrared divergences, also that the magnitude of the quasicondensate 
had to be adjusted such that the jump attained the universal value predicted by Nelson 
and Kosterlitz using renormalization-group methods.! This latter aspect of our work shows 
that the ladder approximation is breaking down close to the critical temperature and needs 
improvement in the critical region. This problem will be addressed in a future publication. 
In this paper we will apply our variational approach developed for the three-dimensional case 
also to the two-dimensional Bose gas, because this approach can in principle also resolve the 
infrared divergences in a manner which is explained in more detail below. 

The best-known experimental realization of a dilute Bose gas in which one can try 
to observe either Bose-Einstein condensation or the Kosterlitz-Thouless transition is spin- 
polarized atomic hydrogen. Although there is no fundamental reason for Bose-Einstein 
condensation not to be attainable in this system, the actual achievement of the critical con- 
ditions is an extremely difficult task due to the inherent instability of the atomic hydrogen 
gas against spin-flip processes which lead to particle loss and in general also cause heating of 
the gas sample. The two-dimensional configuration of course also suffers from these draw- 
backs, but the critical conditions are more favorable than for the three-dimensional case 
and as of the mid 1980’s several experimental groups are therefore pursuing the observation 
of the Kosterlitz-Thouless transition. 0 E In these experiments one uses liquid 4 He films 
to confine the hydrogen atoms in one spatial direction, hence realizing an effectively two 
dimensional Bose gas. 

Other experimental realizations of a dilute Bose gas which are promising candidates for 
the observation of quantum degeneracy effects are laser-cooled alkali atoms and excitons. 
Since evaporative cooling was recently shown to work for magnetically trapped alkalis like Rb 
and A r a,0@ advances towards the observation of Bose-Einstein condensation are particularly 
fast here and recently the actual observation of a Bose-Einstein condensate in an ultracold 
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Rb gas has been reported.ll! Recent developments in MBE-technologies have also led to 
new opportunities to reach the critical conditions for the Kosterlitz-Thouless transition in a 
two-dimensional (dumbbell) excitou gas confined by a double quantum-well structure.^ The 
three-dimensional exciton gas in e.g. pure CU 2 O crystals has a somewhat longer history]^] 
It has even been claimed that the condensation transition has already been observed in such 
a system.il However, mainly due to the lack of a convincing interpretation of the obtained 
experimental data, there is no decisive proof of this fact as yet. 

As mentioned before we will here first study the dilute Bose gas in general by making use 
of a variational approach, and then take spin-polarized atomic hydrogen as a specific exam- 
ple. A variational calculation is a conceptually simple tooSB that has been applied to such 
various subjects as for instance the investigation of the triviality of A 0 4 field theories in differ- 
ent spacetime dimensions, El investigation of the phase diagram of the O(N) vector model, El 
and construction of the ground state of the Sine-Gordon model. El The basic idea is to make 
a Gaussian ansatz for the ground-state wavefunctional and its excitations, which amounts to 
making an ansatz for the dispersion of the quasiparticles in the system, and then to calculate 
the approximate thermodynamic potential = fi tr .+ ('H — 'Htr)tr, where 7 '~t tr is a trial Hamil- 
tonian that has the trial wavefunctionals as eigenstates, Q tr = —ksT Zn(Tr[e _/3( ' WtT ' _AtA ^]) and 
(...)tr is the corresponding grand canonical thermal average. One subsequently minimizes 
this thermodynamic potential with respect to the variational parameters in the Gaussian 
wavefunctionals, and with respect to the expectation value of the held of interest. 

In the case of the dilute Bose gas the results obtained in this manner have to be used in 
the equation of state which makes for example the determination of the critical temperature 
possible. This procedure is a mean-held approach in the sense that no huctuations in 
the expectation value of the held are taken into account. However, it is nonperturbative, 
and clearly goes beyond the one-loop level. The main advantage of using a variational 
approach in two dimensions is that the notorious infrared divergences are automatically 
circumvented because the dispersion relation of the particles will in general display a gap. 
This gap is analogous to the (arbitrary) infrared cutoff k 0 that Popov introduced to define 
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the quasicondensate no(fco)- Note however that in the variational procedure minimization 
of the thermodynamic potential will fix the magnitude of the gap. Hence, there will be no 
adjustable parameter left in the calculation, as was the case in our previous treatment of 
the two-dimensional dilute Bose gas where we had to use the universal jump relation to 
determine the infrared cutoff This convenient feature strongly motivated us to study 
also the two-dimensional system within the framework of a variational calculation. 

Our Gaussian ansatz will introduce a condensate also in two dimensions. Of course, the 
existence of a true condensate in two dimensions is principally forbidden by the Mermin- 
Wagner-Hohenberg theorem.il However, using trial wavefunctionals with a broken symmetry 
in a variational approach can still be usefull, because even if the form of the trial wavefunc- 
tionals deviates considerably from the exact wavefunctionals describing the phase below the 
critical temperature, the estimate for the total energy of the system can be extremely good. 

The organisation of the article is as follows. In section II we derive the desired expression 
for the thermodynamic potential H and we show in particular that minimization with respect 
to the variational parameters automatically introduces the many-body T — matrix into the 
problem. However, it turns out that in its most simple form the variational principle does 
not yield a completely consistent picture. To resolve this problem we have to improve all 
expressions to the T — matrix level. After this improvement we are ready to apply the theory 
to atomic hydrogen. The results are presented in section III. In section IIIA we consider 
the three-dimensional system and section IIIB is devoted to the two-dimensional system. 
Our approach indicates the onset of a new phase in the latter case, which cannot be found 
by an application of perturbation theory around the ideal Bose gas. We are however not 
able to describe the properties of this phase by our variational method but argue that it 
corresponds to the superfluid bound-vortex state. In section IV we draw some conclusions 
from our work. 



5 



II. THERMODYNAMIC POTENTIAL 



We study the dilute Bose gas using the functional integral formulation of the grand 
canonical partition f un ction Thus 



Z gr = e m = J d[ij*}d[iJj]ex.p{-^S[ilj*,ip)} , 



( 1 ) 



where the functional integral is over c-number fields ip*(x,r) and ^(x, r) periodic in imagi- 
nary time over h/3 = h/ksT. The Euclidian action for the dilute Bose gas is given by 

r n/3 



S[VA\ = f Q dr IJ dx -0*(x,r) 



■ d h 2 v 2 

dr 2m ^ 



^(x,r) 



+ \J d *J ^^ x ’ r )^*( x/ ’ r ) V "( x_x, )^( x, ’ r )^( x ’ r )) ’ 



(2) 



with fi the chemical potential and E(x — x') the effectively repulsive interaction potential, 
meaning that the scattering length a is positive. 

Expanding the fields in Fourier modes through 



^(x,r) 



1 

Tnpvyi 2 



k,n 



(3) 



and the complex conjugate expression for ^*(x, r), we can write the action in momentum 
space as 



S[a , a] E( ihuJ n -\-€k d) a k,n a k,n + n y nir E ^ a k+q,n+m a k'-q,n'-rri a k',n , C t k,n • (4) 

k,n ^npV kfc/q 

n,n',m 

In this equation = ti 2 it 2 /2rn and C q = / dx I/(x)e _iq x is the Fourier transform of the 
interaction potential. Since we will be interested in the temperature regime where the spatial 
extension A of the hydrogen atoms is much larger than the range of the interaction potential 
(a/ A <C 1) there is essentially only s-wave scattering taking place in the gas and we can 
neglect the momentum dependence of the interaction potential. Therefore we set E q = Vo- 
The momenta are restricted to either two or three dimensions and V is the area or volume 
of the corresponding system. The bosonic Matsubara frequencies uj n = 27m/h/3 reflect the 
periodicity of the holds in (imaginary) time. 
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We now have the basic ingredients and can start the procedure of variationally calculating 
the thermodynamic potential fh Breaking the symmetry by introducing ^(x, r) = ydTo + 
^//(x, r) and expanding the resulting expression, the action up to quadratic order in first 
instance becomes 

S {2) [a, a*} = -h(3im 0 V + ]rU(3n 2 0 V 0 V + ^ ( -ikuj n + e k - /i + ‘2n 0 V 0 )al n a Kn 

k.n 

1 / 

+ -WqVq^ ( a k,ra a -k,-n T ® k ,n®— k,-n) ) (5) 

" k ,n 

where the prime denotes that k = 0 is excluded from the sum. From this we see that 
the lowest order approximations to the so-called normal and anomalous selfenergies are 
respectively hSn(k, u n ) = 2n 0 Vb and hEi 2 (k, u n ) = no Vo- This gives a strong indication to 
the form of the variational wavefunctional we have to take to correctly describe the dilute 
Bose gas. As the choice of Gaussian wavefunctionals is formally identical to an ansatz for 
the quadratic trial Hamiltonian Ti tr which they are supposed to diagonalize, we now write 
the full action as 



X / 

S[a, a*} = -hPimoV + -h(3nlV 0 V + ^ {-iUu> n + e k - p, + h£n(k,c j n ))a^ n a Kn 

“ k,n 

1 / 

(k, ^ k,— n ^k,ri^— k,— n) 

“ k,n 

+ (2/r 0 Vo — hSn(k, o; n ))^ a k, n a k,n + x( n ok'o — hS 12 (k, cn n ))^ ( a k, n a -k,-n + a k,n<*-k-n) 

k,n k,n 



I n 0 
hf3V 



E 

k',q 

n' ,m 



^O a q,m a k'-q ,n'-m a k',n' + 




£ Coal+q ?n _|_ m ^'q,m^'k,n 

k,q 

n,m 



2 n P v ^ 

n,n f ,m 



^ / O a k +q,n+m a k '-q,n'-m a k , ,n ,a k,n j 



( 6 ) 



where we introduced nontrivial normal and anomalous selfenergies in the first two 
quadratic terms in the right-hand side of Eq. (6) by simply adding and subtract- 
ing Ek, n (^u(k,w„K/ k)n + |hE i2 (k, uj n ) (a k)n od ki _ n + a k , n a_ k _„)). The quantities 
hEu(k, ui n ) and hE 12 (k, u n ) are now considered as the variational parameters of the trial 
wavefunctionals and are determined from a minimization of the thermodynamic potential. 
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As indicated by their arguments, the selfenergies are in principle functions of momentum 
and frequency. However, since we are dealing with a dilute gas in the limit where only s- 
wave scattering is of importance, we will neglect these dependencies and henceforth omit the 
arguments. Note furthermore that we assume the selfenergies to be real, hence hT, i 2 = hE 2 i, 
which amounts to a choice of gauge that is related to our choice of (if;) = yj no being real. 
Performing the Bogoliubov transformation to diagonalize the trial action 

S (tr) [a*, a] = (-^n + e k - // + Mhl)a£ n a k ,n + ( a k,n a -k,-n + °k ,„a_ k -n) , (7) 

k,72 k,n 

we find that the dispersion relation for the Bogoliubov quasiparticles is given by 

nn k = - h') 2 - (hz V2 y , (8) 

where we have introduced fi' — /i — hE n . Moreover, the usual coherence factors!^ for this 
transformation are given by = \[\J (e k — // + hEi 2 )/M\ + yjMl k/ (e k — // + hEi 2 )] and 
t’k = | (e k — yu' + hHi 2 ) / hf2 k — \JhQk/ (e k — /i' + hS 12 )] respectively. The Hugenholtz- 
Pines theorem^ states that the exact normal and anomalous selfenergies satisfy 



/J — hEn(0, 0) — hSi 2 (0, 0) , (9) 

leading to a gapless dispersion. This condition is however not imposed in our variational 
calculation since the Hugenholtz-Pines relation is only valid in a stationary point of the 
actionEi and the variational approach by definition also describes nonstationary points. As 
we will see below, Eq. (9) will nonetheless turn out to hold for the three-dimensional case 
if the temperatures are not too close to absolute zero. In two dimensions Eq. (9) will 
not be valid for temperatures unequal to zero since it would lead to infrared diverging 
integrals in the thermodynamic potential and is therefore excluded from the solution space 
by the variational principle. As mentioned above, this was actually our main motivation for 
applying this approach also to the two-dimensional Bose gas. 

The next step now is to calculate from the functional integral the thermodynamic po- 
tential H which equals 
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( 10 ) 



f! = -j 3 lnZ tr + Ws (i “'>)„. , 

where S^ nt ^ is the interaction term S — S ^ found by a combination of Eq. (6) and Eq. (7). 
Using again the above mentioned Bogoliubov transformation a straightforward calculation 
of this quantity leads to 



, HE 12, n 0 ; /i, T) = ^ngVb - /m + “ e d “ k ) + ^ 77 ^ (^k - £k + //) 



1 

U 



V 



1 1 
2V“ 



k v k 

+ ( 2 n 0 U 0 + //)(rc _ n o) + (^Ei2 - n 0 Vo)hE 12 a 
+ Vo (n — n 0 ) 2 + — — Vo(AtEi 2 ) 2 ar , 



( 11 ) 



where 



i 1 / * \ i 1 V^V ek O' Ar/tn N i e k A 4 V--k 

n = no + — 2 _. ( a k ,n a k,n) = no + — 2 _. ( fc o NyhQk) + 









2MU 



(12) 



k,n k 

is the equation for the density and the dispersion is given by Eq. (8). Furthermore, a 
is introduced as a shorthand notation for 



i i^v,* * i_'i + 2 N(hn k ) 

/ . r'k,ti®-k r n/ ~ y / . 



^12 -' k - U ' k 



2 MU 



(13) 



k ^'<^k 

and iV(a;) = l/(e /3x — 1) is the Bose-Einstein distribution function. 

We now first of all have to minimize this expression for fixed no, chemical potential ji and 
temperature T with respect to the variational parameters AfE 12 and hT,n, or equivalently 
hT .12 and ft! . Putting (dVl/dh'E^) = 0 and (0Q/0 //) = 0 yields respectively 

V’o 2 -nVo + /i' — yU On 



HE 12 = no 



1 + Vo« (1 + Voa)(a + fc£i2 d ^ 12 ) <9AtE 12 



(14) 



and 



dry do/ 

( 2 nl o + // — o)~q^ = [^12(1 + Vo«) — no Vq] , 



(15) 



which are both solved by 



A1E12 = n 0 



Up 

1 + U 0 a 



(16) 
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and 



n' — n — 2nV 0 . (17) 

Calculation of the second derivatives of the thermodynamic potential indeed shows that 
this solution corresponds to a minimum for any Vo > 0. For Vo < 0 the obtained solution 
would correspond to a maximum in the thermodynamic potential, implying an unstable 
solution. This is physically reasonable since a negative Vo implies an effective attraction 
between the constituent particles and thus a negative scattering length a. As a result, a 
magnetically trapped gas indeed turns out to be unstable against the formation of a dense 
(solid or liquid) phase in the (local density) limit that the temperatures are much larger 
than the splittings hu of the one-particle eigenstates in the trap, and the correlation length 
is smaller that the typical trap size.il Quantatively this means that ksT/huj S> 1/ \ a |^> 1, 
where l = \J Ti/moj is the extent of the lowest lying unperturbed harmonic oscillator state of 
the trap. However, in the opposite limit it has recently been shown by a numerical solution 
of the nonlinear Schrodinger equation that a condensate may be metastable for a sufficiently 
small occupation number.il Although this is a very interesting result in the light of the fact 
that a number of alkali atoms, e.g. Cs and Li, have negative scattering lengths, @11 we will 
not consider this highly inhomogeneous case here and therefore restrict ourselves to a > 0. 

Starting with the second equation we immediately conclude that the solution for the 
normal selfenergy is hT n = 2nVo, meaning that the variational calculation yields for this 
particular variable only the one- loop (i.e. Hartree-Fock) expression. The first equation shows 
however, that the anomalous selfenergy hT 12 is obtained in an approximation that goes far 
beyond one loop. In fact, as we will now show, does it constitute the many-body T — matrix 
approximation to the selfcnergy 

The many-body T — matrix T MB (k, k 7 , K; z ), which is diagrammatically depicted in Fig. 
1(a), describes the collision of two atoms having relative momenta k 7 and k before, re- 
spectively after the collision, and a total center-of-mass momentum K and center-of-mass 
energy z in an environment containing the same atoms. The equation for the many-body 
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T— matrix, which by construction sums all the ladder diagrams, can easily be derived if one 
realizes that the free propagation of the atoms is determined by the trial Hamiltonian, that 
can be diagonalized by a Bogoliubov transformation. A somewhat lengthy but straightfor- 
ward calculation gives in the thermodynamic limit 



T MB (k, k\ K; z) = V(k - k') 



dk! 

(2A 



-H(k - k" 



x 



+ 



u\u 2 



v 2 ,v 2 



z — Ml + — z + TlVL i + hVt_ 



U 2 V 2 , 



u 2 + v 2 _ 



z + hQ + — z — hfl . |_ + 

x T MB (k", k', K; z) , 



(1 + N+ + AL) 

(N+-N-)) 



(18) 



where N + = N(hQ + ) and AT_ = The plus sign denotes the momentum argument 

Y + k", and similarly the minus sign denotes the argument ^ — k". This expression for 
the many-body T — matrix has a clear physical interpretation: Two incoming atoms with 
momenta ^ + k" and ^ — k" can scatter out of their momentum states either by the 
creation of two quasi-particles, by the creation of two quasi-holes and by the creation of one 
quasiparticle and one quasi-hole. The net rates for the first two processes are proportional 
to + N + )( 1 + N_) — iV+AA] and v\v 2 _[N + N_ — (1 + N + )( 1 + N_)}, respectively, 

if we take the usual Bose enhancement factors for the scattering into occupied states into 
account. For the last process the quasi-particle can have either momentum y + k" or ^ — k". 
The net rate for these processes are proportional to ■u 2 _Aj_[(l + AL)A+ — AL(1 + N + )] and 
m+^ 2 -[(1 + A^ + )A^_ — Af+(1 + AL)], respectively. Together with the usual energy denominators, 
which are well-known from second-order perturbation theory, this fully explains the structure 
of the right-hand side of Eq. (18). 

Since we neglected all momentum dependencies of the selfcnergies hT,n and hT, u in the 
variational approach, we will also do so in the above equation for the T— matrix. When 
we furthermore set — 1 and = 0 by making use of the fact that the Bogoliubov 
spectrum deviates from — /j , 1 only for a very small momentum interval such that 
e k ~ HT,i 2 < nV o, we find that Eq. (18) reduces to 
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T MB ( 0, 0, 0; z) = V 0 + V 0 T MB ( 0, 0, 0; 2 



dk 1 + 2N(hQ]< 
(27r) d z — 2M2 k 



(19) 



with an error of O(V 0 naA 2 ), which is in general negligible except in an extremely small 
interval of temperatures near absolute zero. By virtue of the fact that we neglected the 
momentum dependence of the scattering potential this equation can easily be solved to give 

Vo 



T MB (0, 0, 0; 0) = 



( 20 ) 



1 + V 0 a 

which finally establishes our claim that the variational calculation gives the anomalous 
selfenergy to within T — matrix approximation, i.e. fr£ 12 = n 0 T AIB ( 0,0, 0;0). Furthermore, 
Eq. (20) implies that a = 1/T MB (0, 0, 0; 0) — 1 /Vq. 

Note that in the calculation we are carrying out, this T — matrix still depends on the 
variational parameters hT,n and fi£ 12 because they enter into the dispersion hfl k . Moreover, 
the many-body T — matrix can be related to the two-body T— matrix, which describes the 
collision of two particles in a vacuum instead of a medium, through^ 

1 1 



T MB (0, 0, 0; 0) T 2B (0, 0; —2//) V 



i N{hn k ) 



Mk 



( 21 ) 



The low-energy limit (e k h 2 /2ma 2 ) of the two-body T — matrix has been studied ex- 
tensively in the literature.^ When configuration space is three-dimensional it is accurately 
represented by 

Airah 2 



T 2B ( 0,0;e k ) = 



m 



( 22 ) 



whereas 



T 2S ((I0;e k ) = 



Anti 2 /m 



(23) 



iri — ln(k 2 a 2 / 8) — 2y 

when configuration space is two-dimensional. In these expressions a denotes the scattering 
length or effective hard sphere radius of the interaction potential and 7 ~ 0.5772 is Euler’s 
constant. These long wavelength expressions will enable us in a later stage to formulate our 
final results solely in terms of the scattering length instead of the full (microscopic) potential 
E(x). 
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Comparing Eqs. (16) and (17) we conclude that the variational calculation as performed 
above leads to an inconsistent theory: The normal selfenergy is obtained to one-loop order 
only, whereas the anomalous selfenergy is obtained in the T— matrix approximation, which 
contains any number of loops. The inconsistency can be better understood by diagrammat- 
ically analyzing the structure of Eqs. (16) and (17), making use of Eq. (12) for the density 
and the definition of a given in Eq. (13). The result is depicted in Fig. 1(b) together with 
the Dyson equation for the normal and anomalous Greens functions. It is clear from these 
diagrammatic equations that iteration will generate the ladder series for the anomalous self- 
energy hT, 12 , but not for the normal selfenergy hT,u. Moreover, it shows that the set of 
diagrams which is being summed in our case of a complex | ijj | 1 theory gives at most the 
T— matrix. This is in contrast to for example the result obtained by Barnes and Ghandour.0 
They show that a similar variational calculation for a real scalar 0 4 theory yields selfenergies 
containing the ladder as well as the bubble diagrams. 

To lift the inconsistency and to obtain an expression for the thermodynamic potential 
independent of the (microscopic) variable Vo we have to improve the derived expressions. We 
proceed as follows. Since it was shown that for any Vo > 0 the solution to the minimization 
procedure yields a minimum in the surface of the thermodynamic potential, we can replace 
it everywhere, except in the condensate contribution \n^Vo for reasons that will become 
clear shortly, by some constant A. Performing again the minimization procedure now leads 
to hE 12 = n 0 A/(l + Aa ) and // = /j — 2nA, where a = 1/T i ' /B (0, 0, 0; 0) — 1/A. It is 
readily seen that taking A equal to the many-body T — matrix evaluated at the minimum of 
12 located at the point in the // — TtL 12 plane given by the above solutions, indeed leads to 
the desired selfcnergies 



HE 12 = n o T MB (0,0,0;0) , 


(24) 


hZ n = 2nT MB (0, 0, 0; 0) 


(25) 



in T — matrix approximation. 

The thermodynamic potential at this minimum in Erst instance equals 
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( 26 ) 



— fi(n 0 ; /.!, T) = -n 2 0 V 0 - fin + -f— l„l I - 



3 -/37iO k \ 



+ 77 77^2 (^k — + aO + (2noT MB (0, 0, 0; 0) + //) (n — ?ro) 

- ' k 

+ T MS (0,0,0;0)(n-n 0 ) 2 , 

from which we still have to eliminate the Vq dependence occurring in the first term, as well 
as the ultraviolet divergence in the fourth term of the right-hand side. Interestingly, these 
two matters are intimately connected, and can both be resolved by writing the sum of these 



two terms as 



o n oVo + - e k + h + 2 ^ k ) - 0(^12) TjE bTf 



V V 2/Lft k 



By also adding and subtracting l(/i£ n) 2 y Jlk and realizing that fiL 12 = 

n 0 T MS (0, 0, 0; 0), we see that Vo is exactly renormalized to the many-body T — matrix if 
one of the T MB ' s in (7iE 12 ) 2 is replaced by Vo, since then 

1<V 0 - \nlV 0 T^(0, 0. 0; 0)^E' 1 + ^ k) = 0. 0. 0; 0) (27) 

due to Eq. (19). The final expression for the thermodynamic potential thus becomes 
ln(n 0 ; li, T) = ^n 2 0 T MB (0, 0. 0; 0) - /in + E HI - 



ily'(RO k - £k + M'+^) + I 



+ V^AE' 



N(hn k ) 



2 V 7 ^ r 2M\ 7 2 V “ 7 V 7 ^ ^k 

+ (2n o T MB (0, 0, 0; 0) + fj!){n - n 0 ) + T MB ( 0, 0, 0; 0 ){n - n 0 ) 2 , 



with fr£ 12 given by Eq. (24) and // = fi — 2nT MB (0 , 0, 0; 0). Note that this derivation by 
no means depends on the validity or violation of the Hugenholtz-Pines relation, and that we 
now have obtained the thermodynamic potential containing only the many-body T — matrix 
and not the interaction parameter Vo. The fourth term on the right-hand side can be shown 
to be at most of 0(Vna 3 ) in three dimensions, and of 0(l//n(l/na 2 )) in two dimensions, 
and will henceforth be neglected. 

Although the above procedure might in first instance seem somewhat heuristic, it has 
indeed led to the correct expression for the thermodynamic potential in T — matrix approxi- 
mation. This we will now show by deriving Eq. (28) directly from the original action S[a*, a] 
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of which we have in Eq. (5) written down the quadratic part. To do so we treat the off- 
diagonal part of S ^ also as a perturbation, and consider only the diagonal part as the 
unperturbed action. Thus (a kn a^ n ) = — Tt/iiTiuJn — e k + (!>')■ In first instance we would then 
immediately write down that 



= ^n 2 0 V 0 - /m 0 + M 1 ~ e /3(ek “ T MB ( 0 , 0 , 0 ; 0)(n - n 0 ) 2 , (29) 

where the last term compensates for the double-counting of the interactions. The first 
term can now be improved to the T — matrix level by performing perturbation the- 
ory in the off-diagonal parts of the action S^ 2 \ The lowest-order term contributes 
— |(ftE 12 ) 2 -p X)k ^ t° the thermodynamic potential, whereas the consecutive terms 

are of O(n 0 a 3 ) resp. O(l//n(l/n 0 a 2 )) and negligible. Combining this lowest-order term with 
\uqVo the thermodynamic potential now reads 

iff = \n 2 0 T MB { 0 , 0 , 0 ; 0) - /m 0 + ^X>( 1 - - T MB { 0 , 0 , 0 ; 0)(n - n 0 ) 2 . (30) 

This result now has to be compared with Eq. (28). Firstly we have to show, lest they be 
equal, that 



ksT vV '-'’ r wn ") + 



V 



-J2 H 1 



_ k B T' ’ _ -/3(e k -/ i ')'j 

hV nn k v ^ [ J 



(31) 



k ^ ^ k '"“k v k 

That this equality holds to a very good approximation can be demonstrated simply by 
performing an expansion of hfl kl which leads to 



ln{ 1 - e -^ k ) = ln( 1 - e^”')) - (5 



{KZ 



12 



2(e k - 



-A^(e k - //) + h.o.t. 



(32) 



where the ratio of the higher order terms and the Erst order term in the expansion is either 
of 0(noaA 2 ) or of 0((noA 2 /ln(l/noa 2 ))) in the regime where these parameters are small. 
Close to zero temperature these parameters become large, but what has to be compared 
in that situation is the contribution of these terms to the total thermodynamic potential. 
Due to the Bose factors their contribution becomes negligible. Secondly, the two respective 
terms (n — no) 2 in Eqs. (28) and (30) also differ clue to the difference in dispersion. Again 
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expanding Mik, one can show that the error here is of the same order as above and negligible. 
Thus Eq. (28) obtained from the modified variational calculation is indeed the correct result 
for the thermodynamic potential. 

The above argument shows that the results obtained in Eqs. (24), (25) and (28) can 
also be interpreted in a different manner. First we break the symmetry of the theory by 
hand. We subsequently calculate the normal and anomalous selfenergies hT,n and hT, u 
in the T— matrix approximation. This straightforwardly gives Eqs. (24) and (25) (see also 
Ref. 12). We then also calculate the thermodynamic potential in T— matrix approximation, 
yielding as we have seen Eq. (28). Thus we again arrive without reference to a variational 
calculation at the same set of equations and in particular at the same thermodynamic 
potential as a function of condensate density, chemical potential and temperature. 

The final step we now have to take is to find the value of n 0 which minimizes the 
thermodynamic potential at given p and T. After this the properties of the system are 
fully determined and it is possible to calculate e.g. the condensate fraction as a function 
of temperature and the superfluid density. Using the fundamental thermodynamic identity 
E = fl + TS + pN, with the entropy S given by 

s = ~k B Y,i N(h{i k )in N (hn k ) - (i + N(nn k ))in ( i + N(nn k ))] , (33) 

k 

we can also calculate the specific heat from Cy = (|^ ) N ■ 

Because this program has to be implemented numerically, we now restrict ourselves to 
the specific example of spin-polarized hydrogen atoms in both a three-dimensional (section 
11IA) and a two-dimensional (section I11B) configuration. In the case of the system being 
three-dimensional, it leads to improvements on existing approaches to this problem, and we 
therefore in particular present results for the condensate fraction, the specific heat and the 
p — n _1 diagram. 
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III. SPIN-POLARIZED ATOMIC HYDROGEN 



In order to obtain Bose-Einstein condensation one magnetically traps hydrogen atoms 
in the doubly polarized | d) =|T^) state of the upper hyperhne manifold, resulting in a 
well-known effectively repulsive triplet interaction between the atoms with a value of the 
scattering length of 1.34ao-S Typical sample densities are 10 13 to 10 14 atoms per cubic 
centimeter with a corresponding critical temperature on the order of tens of microkelvins. 

For the two-dimensional system one uses atoms in the also doubly polarized j b ) =||-^) 
state of the lower hyperhne manifoldi adsorbed on a superfluid helium him. The effective in- 
teraction between the adsorbed atoms is in a very good approximation given by the weigthed 
average of the volume triplet potential, using the square of the bound state wavefunction 
pertaining to the motion perpendicular to the surface, as a weighting factor.ll This pro- 
cedure leads to a two-dimensional scattering length with a value of 2.40a 0 .il The densities 
one aims at in this experimental setup are typically some 10 13 atoms per square centimeter, 
which corresponds to a critical temperature in the range of 50 to 100 millikelvin. 



A. The three dimensional case 

In Fig. 2 we show for three different values of noa 3 the curves in the hT 12 — p 7 plane 
on which the equality HT 12 = noT MB (0, 0 , 0 ; 0) holds, as well as two curves on which 
the chemical potential p is constant, both at a temperature of 30 fiK. The straight line 
fjf = —hT 12 corresponds to the conditions where the Hugenholtz-Pines theorem is satisfied. 
The area T1T12 >\ n' | is of course not part of the solution space since it would imply a 
complex (unstable) dispersion relation. By calculating the thermodynamic potential at the 
intersections of the lines we infer the relation between and the condensate density no for 
a given chemical potential. 

This is shown in Fig. 3 for the same two values of the chemical potential. From this plot 
we find that the no minimizing the thermodynamic potential for given /i indeed corresponds 
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to the situation where the Hugenholtz-Pines theorem is satisfied, i.e. // — Note 

that this is a result of the minimization procedure and is not assumed beforehand. The 
Hngenholtz-Pines relation however turns out not to hold in every situation. For extreme 
conditions, very close to zero temperature, the position of the minimum is shifted slightly 
from the Hugenholtz-Pines line. (See also Fig. 13 in the next subsection.) At T = 0 one 
can, starting from Eq. (28) with general hT, 12 and //, prove also analytically that this is the 
case. The deviation from the Hugenholtz-Pines condition is due to an overestimate of the 
effective interaction between the atoms. Indeed, a smaller value of the T— matrix, together 
with a larger value of the condensate density to ensure that the chemical potential /1 remains 
the same, shifts the minimum towards the Hugenholtz-Pines line. This overestimate of the 
effective interaction is due to the neglect of the influence of the coherence factors in the 
expression for the many-body T— matrix in Eq. (19), i.e. due to putting w k = 1 and u k = 0. 
Taking these coherence factors into account will cause a decrease of the T— matrix which 
presumably solves the problem that the minimum of the hi — no— curve does not coincide 
with the Hugenholtz-Pines criterion. The region where the neglect of the coherence factors 
becomes important corresponds to the condition naA 2 ~ 1, which is indeed in the extremely 
low temperature range where our numerical results show a minimum in the thermodynamic 
potential that is slightly shifted from the Hugenholtz-Pines line. Since this temperature 
range will be very difficult to obtain experimentally, we will not consider it further here. 

Outside this range, all the properties of the system can now be calculated accurately. In 
Fig. 4 we plot the condensate fraction as a function of temperature for a density of 10 14 cm~ 3 . 
The inset shows the behavior near the critical temperature following from our calculation, 
as well as from a calculation in which we have replaced the many-body T — matrix by the 
two-body T— matrix. The latter is essentially what has been done in earlier approaches to 
the dilute Bose gasJ Using the two-body T— matrix however produces a double-valuedness 
in this curve, leading inevitably to the conclusion that there is a first-order phase transition. 
This double-valuedness also shows up in the equation of state depicted in Fig. 5. Again, 
the equation of state is single-valued if the many-body T— matrix is used. Our variational 
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approach has therefore succesfully resolved the long-standing problem concerning the double- 
valuedness in the equation of state. It shows that it is essential to include the effect of the 
surrounding medium on the collision of two hydrogen atoms, which leads effectively to a 
temperature dependent scattering amplitude and a softening of the dispersion near the 
critical temperature. We calculated the effect of the medium on the effective scattering 
length explicitly for the above situation where n = 10 14 cm -3 , and show in Fig. 6 the 
ratio of the many-body T — matrix to the two-body T — matrix as a function of temperature. 
The fact that the many-body T— matrix, and therefore the effective interaction of the atoms, 
goes to zero at the critical temperature is correct when neglecting momentum and frequency 
dependencies as can be shown from an exact renormalization-group equation for the four- 
point vertex function.® 

ffowever, the many-body T— matrix T MB (k, k', K; z) does not go to zero for nonzero 
momenta and frequencies. This explains the feature seen in the inset of Fig. 7 , where we plot 
the p—n~ l diagram resulting from our calculations. Near the critical density the curve shows 
a van der Waals loop, which, by means of a Maxwell construction, would still imply a first- 
order phase transition. This behavior is due to the fact that the thermodynamic potential 
Eq. (28), and hence the pressure p = —Q/V, always contains the many-body T — matrix 
evaluated at zero momenta and energies. Therefore, the van der Waals loop shown in Fig. 
7 is solely a consequence of the behavior of the ratio T A/S (0, 0, 0; 0)/T 2S (0, 0 ; 0 ) near the 
critical temperature (see Fig. 6). It would not have been present had we used a more accurate 
expression for the thermodynamic potential, which in principle contains integrals over the 
many-body T — matrix evaluated at finite momenta and energies. Hence, the first-order 
phase transition implicate from Fig. 7 is an artefact of our long wavelength approximation 
which breaks down near the critical temperature. We stress however that the softening of 
the dispersion discussed earlier, resolving the problems related to the double-valuedness in 
the equation of state near the critical temperature in contrast is not an artefact because the 
dispersion relation does contain T AIB ( 0, 0, 0; 0 ) in the long wavelength limit. 

Finally in Fig. 8 we show the specific heat Cy as a function of temperature, again for 
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a density of 10 14 cm~ 3 . Close to T c our results in principle show an artificial discontinuity 
which is due to the first-order phase transition we find as explained above, and is therefore 
not explicitly shown. This concludes our discussion of the three-dimensional system, and 
we now turn to the two-dimensional case. 

B. The two dimensional case 

In Fig. 9 we show the solution curves of = n 0 T MB (0, 0 , 0 ; 0 ) for three different 
values of rioa 2 at a temperature of 50 mK, as well as two curves on which the chemical 
potential (i is constant. In addition to the area hT, i 2 >| /i' |, also the Hugenholtz-Pines line 
itself is not part of the solution space since it would imply a gapless dispersion and therefore 
infrared divergences. This is also the cause of the essential difference of Figs. 9 and 2. To 
avoid infrared divergences in the first case, the solution curves of 7i£ 12 = n 0 T MB ( 0 , 0 , 0 ; 0 ) 
are forced to start in the origin of the // — frS 12 plane, in contrast to what happens in three 
dimensions (cf. Fig. 2). This is a crucial point of difference and results in the fact that the 
curves of the thermodynamic potential versus chemical potential all start in one and the 
same point, corresponding to /i = ft£ 12 — fjf — 0, as shown in Fig. 10, where n 0 increases 
going from curve 1 to curve 4. 

From this graph the following is then immediately clear: For < 0 the value of no yielding 
the lowest thermodynamic potential is zero, whereas for fi > 0 the value of no has to be 
infinite to give the lowest value of Thus for yU > 0 we cannot fix the number of particles in 
the system to a finite value, and we conclude that the variational approach signals the onset 
of an instability at yU = 0, but cannot describe the phase below the critical temperature. 
Above T c it describes the two dimensional Bose gas in the ladder approximation. What we 
can do is to deduce from the condition /i = 0 the relation between the critical temperature 
and density as it follows from our variational approach. This is shown in Fig. 11, where we 
plot n c versus T c . As expected we find that the critical temperature is overestimated, and 
n c JS? c is smaller than 4. This is consistent with the fact that the variational approach is of a 
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mean-field type, since no fluctuations in the background field have been taken into account. 
Fluctuations, and in particular phase fluctuations associated with vortices, should reduce 
the real critical temperature and result in the universal jump of the superfluid density. This 
is schematically shown in Fig. 12, where we plot the behavior of the superfluid density versus 
temperature. Only for temperatures larger than T C MF , which is itself larger than T BT at 
which the Nelson-Kosterlitz criterion n s A 2 = 4 is satisfied, can we apply our theory. 

We would like to point out that T C MF cannot be found by a straightforward application 
of perturbation theory around the ideal Bose gas. Indeed, calculating the free energy F and 
the (normal) selfenergy in the ladder approximation would result in 

// = n - hTiu = n- 2nT MB (0, 0 , 0 ; 0) (34) 

and in the equation of state 

n = J2 N ( e k “ tf) » ( 35 ) 

k 

where the latter is now obtained by a minimization of the free energy. At a fixed temperature 
and density, we first of all determine //' from Eq. (35) and subsequently // from Eq. (34). 
Clearly, this calculation can be performed at any temperature and we just obtain the ideal 
gas result that a phase transition occurs only at zero temperature. Notice however that for 
temperatures larger than T C MF the above procedure of minimizing the free energy F leads, as 
expected, to exactly the same results as our minimization of the thermodynamic potential. 

Finally we show in Fig. 13 the thermodynamic potential as a function of condensate 
density for a fixed chemical potential in the case of zero temperature where a real con- 
densate can exist in the two-dimensional system. In this case the solution curves of 
hEi 2 = noT AIB (0, 0 , 0 ; 0) do not start out in the origin of the /i' — hSi 2 — plane, but es- 
sentially behave as depicted in Fig. 2, because without thermal fluctuations there are no 
infrared divergences when jJ = — hE 12 . Again in this extreme situation we find that the 
Hugenholtz-Pines relation is not satisfied, which, as in three dimensions, is due to the ne- 
glect of the coherence factors u k and Vk in the equation for the many-body T— matrix by 
our variational approach. 
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IV. CONCLUSIONS 



We performed a variational calculation for a weakly-interacting Bose gas and applied 
this method to spin-polarized atomic hydrogen in both the three-dimensional and the two- 
dimensional case. We showed that the variational approach does not automatically lead to 
a consistent description of the system under study, and that the number of diagrams that is 
being summed can be rather restricted. In particular, the normal selfenergy is obtained in 
one-loop approximation only. The anomalous selfenergy contains the ladder diagrams. We 
consequently improved all relevant expressions to the T — matrix level. 

In three dimensions this approach resolves the problem of double- valuedness in the equa- 
tion of state, that leads to the conclusion of a first-order phase transition and is found in all 
previous approaches to the dilute Bose gas. The essential ingredient to resolve this problem 
is the use of the many-body T— matrix, which takes into account the fact that two particles 
collide in a medium, and not in a vacuum. The (small) van der Waals loop we still find 
in the p — n~ l diagram is a result of the neglect of the momentum and energy dependence 
of the T— matrix. The simplification in the coherence factors u k and Uk leads to an over- 
estimate of the effective interaction, and explains the fact that in extreme conditions, very 
far from the critical temperature, the Hugenholtz-Pines relation is not satisfied. Using the 
full many-body T — matrix T AIB (k, k', K; z) will resolve these issues. Not too close to the 
critical temperature or to zero temperature the theory gives an excellent description of the 
dilute Bose gas. We also calculated condensate fraction and specific heat as a function of 
temperature. 

In two dimensions the variational approach, having a mean-field nature, signals only the 
onset of an instability in the theory, but unfortunately cannot describe the phase below the 
critical temperature. This is connected to the fact that within this variational approach it 
is not possible to find solutions where the dispersion relation of the quasiparticles displays 
a gap, playing the same role as the infrared cutoff in the theory with a quasi- condensate 
and enabling us to describe the system also below the critical temperature. The critical 
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temperature we can infer from the occurence of the instability is an overestimate of the true 
critical temperature due to the neglect of the crucial phase fluctuations. In this respect it is 
important to note that introducing a quasi-condensate no(ko), followed by a diagrammatic 
calculation of the equation of state, 0 does give correct information about the Kosterlitz- 
Thouless phase and thus incorporates the essential physics of the superfluid state. This in 
contrast to the mean-held variational approach presented here. 

Finally we want to remark that solving the problem of the dilute Bose gas using the 
full T — matrix Eq. (18) is a highly nontrivial task due to infrared difficulties even in three 
dimensions. A renormalization-group study however will solve these difficulties and more 
importantly will also correctly describe the critical region. A study along these lines is 
underway. 
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FIGURE CAPTIONS 



Fig. 1 Diagrammatic representations. 

(a) Diagrammatic representation of the T— matrix equation. 

(b) The normal and anomalous selfenergies as found in the variational approach together 
with the expressions for the exact Green’s functions. Note that iteration of the equation for 
E 12 generates the ladder diagrams. This is not the case for Sn. 

Fig. 2 Lines of constant chemical potential ((4) and (5), the first having the largest chemical 
potential), and solution curves of the equation hEi 2 = noT AIB (0, 0 , 0 ; 0) in the //' — hEi 2 
plane for (1) n 0 a 3 = 4.8 10~ 12 , (2) n 0 a 3 = 5.0 10 -12 and (3) n 0 a 3 = 5.2 10~ 12 . Also shown 
is the line on which the Hugenholtz-Pines relation is satisfied. 

Fig. 3 Thermodynamic potential versus condensate density for the same two values of the 
chemical potential as in Fig. 2 and at a temperature of 30/iK . The dots indicate the starting 
points on the Hugenholtz-Pines line of the corresponding curves in the /. l ' — hE 12 plane. The 
graph shows that at the minimum the Hugenholtz-Pines condition is indeed satisfied. 

Fig. 4 The condensate fraction as a function of temperature for a total density of 10 14 cm~ 3 . 
The inset shows the behavior of the condensate fraction near the critical temperature reveal- 
ing the double-valuedness (curve (1)) if one uses the two-body T— matrix, which is resolved 
if one uses the many-body T— matrix that includes the effect of the medium on the collision 
(curve (2)). 

Fig. 5 The equation of state for small // at the critical temperature for a system with 
n = 10 14 cm -3 , showing the double-valuedness discussed also in Fig. 4. 

Fig. 6 The ratio of the many-body T— matrix to the two-body T— matrix at a density of 
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10 u cm 3 , showing that T MB (0 , 0 , 0 ; 0) goes to zero at the critical temperature. 

Fig. 7 Pressure versus inverse density at T = 50 fiK. The inset shows a small van der 
Waals loop present around the critical density (see text). 

Fig. 8 Heat capacity versus temperature when n = 10 14 cm“ 3 . Note that 3n/cs/2 corre- 
sponds to 0.002071 

Fig. 9 Lines of constant chemical potential ((4) and (5), the first having the largest chemical 
potential), and solution curves of the equation hS 12 = n o T AIB (0, 0 , 0 ; 0) in the //' — h£ 12 
plane for (1) n 0 a 2 = 10” 5 , (2) n 0 a 2 = 5 10 -5 and (3) n 0 a 2 = 10 -4 . Also indicated is the 
Hugenholtz-Pines line. 

Fig. 10 Thermodynamic potential versus chemical potential at a temperature of 50 mK and 
(1) n 0 a 2 = 0, (2) n 0 a 2 = 1CT 5 , (3) n 0 a 2 = 10~ 4 and (4) n 0 a 2 = 5 10” 4 . The dot indicates 
the starting point of the corresponding curves in the fi' — hS 12 plane, which lies in the origin 
for all values of noa 2 . 

Fig. 11 Critical density versus temperature found from the criterion fi = 0, the point at 
which the variational approach signals the occurrence of an instability. 

Fig. 12 Schematic behavior of superfluid density versus temperature. The full line results 
from an RG-analysis, that leads to the exact critical temperature T ^ T , obeying n s A 2 = 4. 
The dashed curve corresponds to a mean-held result with a critical temperature T A,r . Below 
T aif the variational approach does not give information obout the properties of the system. 

Fig. 13 Thermodynamic potential versus condensate density for fixed chemical potential 
at zero temperature, showing that the Hugenholtz-Pines condition, valid only at the end- 




point of the curve, is not satisfied in the minimum of this curve. 
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Figure 1(b) 






